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ABSTRACT 



Measurements of a continuous hydromagnetic micro- 
pulsation event, which occurred in the pre-dawn 
magnetosphere on July 29, 1979, are reported as observed by 
wave and particle instruments onboard the P78-2 (SCATHA) 
spacecraft. Calculations of the Poynting vector make it 
clear that the wave is a traveling Alfven wave guided by the 
earth's magnetic field line. Plasma densities are 
calculated at L = 7. The phase relationship between the 
plasma flux and the electromagnetic fields is investigated. 
A self-consistent set of particle density and particle 
velocity data is determined, in order to verify the 
measurement of the electric field. Current magnetospheric 
models predicting Alfven wave speeds, Alfven wave periods 
and plasma densities are supported by the study. 
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I. INTRODUCTION 



A. GEOMAGNETIC MICROPULSATIONS 

Geomagnetic micropulsations are low frequency 
perturbations of the Earth's magnetic field. They are 
generally short-lived variations whose amplitude is small 
.when compared to the magnitude of the Earth's main field. 
In contrast to long term variations of the Earth's magnetic 
field which are of internal origin, geomagnetic 
micropulsations, as well as magnetic storms, are thought to 
be generated by solar forces. [Ref. 1] 

The gamma, defined as. 10“^ Tesla, is commonly used when 
measuring these extremely small magnetic fields. 
Micropulsations have amplitudes ranging from fractions of 
gammas to tens of gammas and periods ranging from 
approximately 0.1 second to 10 minutes. Micropulsations may 
be divided into two main classes. Those with regular or 
continuous patterns are called Pc's and those with irregular 
patterns Pi's. The system used in this paper for 
categorizing continuous pulsations is described in Table 1. 
[Ref. 1] 

B. CONTINUOUS PCS OSCILLATIONS 

Micropulsations with the longest periods, PcS's, may 
have amplitudes near 100 gamma, but generally are seen to 
have amplitudes of the order of ten to twenty gamma. A 
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group of Pc5 ' s may last from 10 minutes to several hours. 
Pc5 ' s have been noted to be absent on days of very low 
magnetic activity and on very disturbed days. They 

generally occur during the recovery periods of several days 
which follow moderate magnetic disturbances. [Ref. 1] 

TABLE 1 

CONTINUOUS PULSATION CLASSIFICATION 
Type Period ( sec . ) 



Pci 


0.2 


- S 


Pc2 


S 


- 10 


Pc3 


10 


- 4S 


PC4 


4S 


- ISO 


PCS 


ISO 


- 600 



1 . Previous Observations 

Magnetic field micropulsations have been observed 
from ground stations since the 19th Century. In 1861 
Stewart observed what he described as a great magnetic 
disturbance at Kew Observatory [Ref. 2]. Since then 
thousands of papers have been published concerning 
micropulsations. Ground station characteristics of PcS's 
were summarized by Barfield and McPherron [Ref. 3]. Ground 
stations report PcS's most frequently at dawn and dusk. 
These data indicate that it is common for PCS ' s to span 60 
degrees of longitude. Simultaneous observations by magnetic 
conjugate ground stations suggested a generation mechanism 
located near the equator. 

With the advent of satellite exploration in the late 
1960 's the understanding of long period micropulsations was 
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furthered with data taken by satellite magnetometers. 
Several studies were made of data taken by spacecraft in 
geosynchronous orbit. For example, correlations were 
established between micropulsation observations in 1967 from 
Applied Technology Satellite 1 (ATS-1) and ground station 
data [Refs. 3 and 4]. This analysis showed that the onset 
of 85 percent of the events occurred during the expansion 
phase of magnetospheric substorms, suggesting that Pc5 ' s 
could be generated by storm enhanced ring currents . 

In the early 1970 's plasma scientists began to 
correlate data made available by multiple sensors on 
spacecraft. Comparison of magnetic field data and high 
energy plasma data, from Explorer 45, with ground station 
data showed that Pc5 resonance regions existed in the 
magnetosphere at high altitudes [Ref. 5]. It was suggested 
that these resonance regions could be perturbed by a drift 
mirror instability. A good review of plasma instabilities 
is given by Chen [Ref. 6]. 

ATS-6 provided an early example of combined magnetic 
field and plasma flow measurements in 1974 [Ref. 7]. Using 
multiple sensors it was possible to calculate the Poynting 
vector using measured magnetic fields and inferred electric 
fields. This case study showed the existence of a standing 
Alfven wave. 

Multi-sensor observations have become more common in 
the past decade. The addition of mass spectrometers brought 
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new dimension to magnetospheric research. A recent example 
is a Pc5 observed by Dynamics Explorer I in July 1982 [Ref. 
8]. Using a mass spectrometer, electric field and magnetic 
field measurements, it was possible to compare simultaneous 
measurements of the E X B drift velocities and particle 
derived plasma drift velocities during a Pc5 event. The 
plasma was shown to be rotating perpendicular to the 
magnetic field. Complete analysis of the plasma data showed 
that significant quantities of singly ionized helium, oxygen 
and nitrogen, as well as hydrogen, were present in the 
plasma, which was important in matching plasma theory and 
observations . 

GEOS-2 data from the double probe electric field 
experiment have been reported for PcS's observed during 1979 
[Ref. 9]. Observed correlations between wave amplitude and 
plasma density data were used to support current hydro- 
magnetic wave period models. 

Multi-satellite observations can be used to address 
the spatial structure of PcS's. Combined magnetic field 
data from the SCATHA, GOES 2, GOES 3 and GEOS 2 satellites 
were used to examine the structure of a field aligned Pc5 
observed in November 1979 [Ref. 10]. The study showed that 
the wave had an antisymmetric standing structure along the 
field line, and that the wavelength was much shorter than 
the field line. These are characteristics that have not yet 
been explained by existing theories . 
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Analysis of these observations, leading to the 
determination of wave modes, growth characteristics and 
propagation characteristics, enhance the understanding of 
basic plasma physics and those processes which determine the 
dynamics and structure of the magnetosphere. 

2 . Theory 

H. Alfven showed, in 1942, that low frequency 
electromagnetic waves can propagate in fluids of high 
conductivity in the presence of a constant magnetic field 
[Ref. 11]. 

These waves are difficult to observe on a laboratory 
scale. S. Lunquist determined that for a strong interaction 
to occur between electromagnetic and hydrodynamic systems 
the criterion which must be met is 

BLo(y/p)l/2 >> 1 (1.1) 

where B is the magnetic field, L is the dimension of the 
conducting fluid, o is the conductivity and p is the density 
of the fluid [Ref. 12]. Since the dimensions of the 
magnetosphere are large, this type of interaction is 
possible . 

In a Newtonian fluid, of perfect conductivity and 
density p, small perturbations of the magnetic field, Bq, 
will propagate along the magnetic field line with velocity 
given by 

Va = co/k = Bq / (yoP)l/2 (1.2) 
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which is called the Alfven velocity [Ref. 6]. 

In 1952, Alfven compared hydromagnetic waves to the 
vibrations of an elastic string with its ends fixed at 
magnetic field conjugate points [Ref. 13]. A good summary 
of the analogy is given by Sugiura and Wilson [Ref. 14]. 
The derivation of this analogy is described showing that the 
plasma acts as the mass of the string. The calculation of 
perturbation modes of these string vibrations, as applied to 
the earth's magnetic field, is presented. 

Hasegawa and Chen separate the theory of long period 
micropulsations into two categories and set forth a theory 
for each type [Ref. 15]. One type can be explained as a 
resonant Alfven wave excited at a local field line by a 
Kelvin-Helmholtz instability at the magnetopause. The 
second type is interpreted as an excitation of a surface 
eigenmode, a Kinetic Alfven Wave, excited by an impulse 
energized ring current. 

In this paper a Pc5 event is described which 
exhibits some characteristics of a classical Alfven wave. 
Its source is believed to be the ring current which 
generates a traveling wave along the magnetic field line. 

C. COORDINATE SYSTEMS 

The coordinate systems pertinent to this analysis of the 
Pc5 data are shown in Figure 1. The magnetic field com- 
ponents (Bx,By,B 2 ) are described in topographic coordinates. 
X is defined as North, Y points East, and Z points down. 
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SOLAR ECLIPTIC 
SYSTEM 




Figure 1. Coordinate Systems 
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toward the center of the earth, in the radial direction. 

The electric field data are despun into two components 
(Ejj and E 2 ) of the solar ecliptic coordinate system. The X 
direction in this system is parallel to the earth-sun line 
and points in the sunward direction. The Y axis is perpen- 
dicular to the earth-sun line and points toward dusk. The Z 
axis is right-hand perpendicular to the X-Y plane. 

The satellite coordinate system is a rotating coordinate 
system with X pointing in the direction of the satellite 
spin axis, which is essentially held constant throughout the 
orbit. The spin axis is nearly anti-parallel to Y (SE) . The 
Y and Z axis rotate through the plane perpendicular to the 
spin axis. 

D . SCATHA SATELLITE 

The U. S. Air Force P78-2 satellite was launched into a 
near-geosynchronous orbit in January 1979 as part of a joint 
NASA/ Department of Defense program to study Spacecraft 
Charging At High Altitudes (SCATHA). A schematic of the 
SCATHA vehicle is shown in Figure 2. 

1 . SCATHA Orbit Parameters 

The SCATHA orbit is elliptical with a perigee of 5.5 
Re and an apogee of 7.7 Re. The inclination of the orbit is 
7.5 degrees and the period is 23.5 hours. The satellite 
spin period is 59.5 seconds, with a spin axis in the orbit 
plane and perpendicular to the earth- sun line. [Refs. 16 and 
17] 
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UCSO PUVSMA DETECTOR 




SCATHA SPACE VEHIClf (P78-2) 



Figure 2. Schematic of the SCATHA Satellite 
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Figure 3 shows the satellite's radial distance from 
the earth's center in earth radii, (1 Re = 6378 km), for the 
24 hours of Day 210. The Pc5 is observed from 4 to 7 hours 
UT during a continuous decline in SCATHA's radial distance. 
Figures 4 and 5 show the geographic latitude and longitude 
during the same period. During the Pc5 event the geographic 
latitude and longitude remain nearly constant. The 
variation in SCATHA's latitude is a result of the 7.5” 
inclination in the orbit. The variation in longitude 
reflects the fact that the orbit is not geosynchronous, but 
elliptical. The result of the slighly inclined elliptical 
orbit is a continuous change in the satellite's geographic 
longitude as well as it's velocity and radial distance from 
the earth. SCATHA's velocity varies from 2.55 km/s, at 
apogee, to 3.74 km/s, at perigee. The virtue of this 
slightly elliptical orbit is that the satellite can sweep a 
wide sector of the earth's magnetosphere while introducing 
only small errors in the raw magnetospheric data. 

In the analysis of magnetospheric data, the satel- 
lite's location relative to the earth's magnetic field is 
important in the classification of observed phenomena. A 
first approximation of the shape of the earth's magnetic 
field is that of a dipole given by 

R = r sin^e (1.3) 
where r is the distance from the earth's center to the 
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RADIAL DISTANCE 



SCATHA DAY 210 




Figure 3. Radial Location of SCATHA on July 29, 1979 
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GEOGRAPHIC LATITUDE ( DEGREES NORTH 



SCATHA LATITUDE 




Figure 4. Geographic Latitude of SCATHA on July 29, 1979 
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GEOGRAPHIC LONGITUDE ( DEGREES EAST 



SCATHA LONGITUDE 




Figvire 5. Geographic Longitude of SCATHA on July 29 , 1979 
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intersection of the dipole field line with the equatorial 
plane. R and 9 are the polar coordinates of the magnetic 
field line. The parameter used to define the location of a 
magnetic field line is the Mcllwain L [Ref. 18]. Using the 
dipole approximation for the magnetic field, the L of the 
field line is equal to the number of earth radii at which 
the field line crosses the equatorial plane. The value of L 
at other latitudes is given by 

L = Re / cos^Gj^ (1*4) 

where is the magnetic latitude. 

A view of the magnetic field dipole approximation is 
shown in Figure 6 where field lines for L = 5, 6, 7, and 8 
are plotted. Scatha's dipole magnetic position between 0400 
UT and 0700 UT is plotted. At 0400 UT SCATHA was at L = 
8.03 and descended to L = 6.59 at 0700 UT. It is during 
this period that the Pc5 was observed . 

The dipole approximation for the earth's magnetic 
field is most exact at dawn and dusk. The approximation is 
fairly inexact at noon (where the field is compressed) and 
poor at midnight (where the field is extremely elongated). 
The distortion of the earth's magnetic field is largely due 
to the solar wind. A more exact description of the shape of 
the magnetic field is given by Hess [Ref. 19]. 
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SCATHA POSITION AT 0700 UT 



O SCATHA POSITION AT 0-^00 UT 
Figure 6. Mcllwain L and Geomagnetic Latitude of the Pc5 
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2 . Instruments 



The spacecraft is a cylinder approximately 1.75 
meters in length and diameter. The instruments used in this 
study are the NASA Goddard Space Flight Center electric 
field monitor, the NASA Goddard Space Flight Center mag- 
netometer and the University of California, San Diego (UCSD) 
charged particle experiment. The locations of these 
instruments are shown in Figure 2. [Ref. 17] 

a. The NASA Goddard electric field detector 
consists of two antennas that are 50 meters in length 
projecting into the spin plane which are used as floating 
probes. The antenna wire is composed of beryllium copper. 
Each antenna is insulated except for the last 20 meters, 
which is the active part. Differential signals between the 
antennas give the ambient electric fields . The minimum 
measurable field is 0.05 mV/m. The time resolution of the 
electric field measurements is limited by the satellite spin 
period (59.5 seconds). 

SCATHA data are despun into the X and Z 
components of the solar ecliptic (SE) coordinate system. By 
this method the two components of the measured electric 
field (E^ and E^) , which are in the plane perpendicular to 
the satellite spin axis are known. It will be shown that 
the component of the field data reflects the component of 
the Pc5 electric field which is perpendicular to the earth's 
magnetic field when the satellite is near local dawn. The 
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transformation from the measured fields to the ambient 
fields requires the addition of a 0.3 mV/m sunward field to 
account for satellite motion. [Refs. 16 and 17] 

b. The GSFC magnetometer is a triaxial flux gate 
magnetometer mounted on a 4 meter boom. The resolution is 
0.3 gamma with four vector measurements per second. Data 
from the GSFC magnetometer are used in raw form and despun 
into the three vectors listed of the topographic coordinate 
system. [Refs. 16 and 17] 

c. The UCSD charged particle experiment is an 
electrostatic analyzer that measures both ion and electron 
flux. Particle data from the electrostatic analyzers (ESAs) 
are used to measure the ion and electron distribution 
functions and deduce the particle density, temperature and 
velocity. The UCSD particle detectors consist of five ESAs 
arranged in two rotating pairs of electron and ion 
detectors, and one fixed ion detector. All three of the 
sensors of this instrument are contained in a single package 
mounted on the forward end of the satellite. [Refs. 16 and 
17] 

One of the pairs of rotating detectors covers 
the energy range from 1 eV to 81 keV (the N/S detectors). 
The other pair (the E/W detectors) and the fixed detector 
cover the 1 eV to 1800 eV energy range. The analyzers have 
energy resolution of 20% (6E/E) and require 16 seconds for 
one 64 step energy scan. [Refs. 16 and 17] A more complete 
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description of the satellite and instrumentation is 
contained in the IMS Source Book [Ref. 20] and the P78-2 
Technical Report [Ref. 21]. 
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II . OBSERVATIONS 



A. PARTICLE AND FIELD DATA 

During the morning hours of Day 210 SCATHA recorded a 
disturbance in the earth's magnetic field, beginning at L = 
7.96, as the spacecraft approached the maximum inclination 
of its orbit, 7.5 degrees. A strong magnetic field 
micropulsation, a Pc5, was recorded on each of SCATHA 's 
sensors . 

Figure 7 is an equatorial plane plot of SCATHA 's orbit 
on July 29, 1979 (Day 210). The location of the extremes of 
the SCATHA orbit are shown with respect to the earth sun 
line. The sector of the orbit, in which the Pc5 event was 
observed, is also indicated in Figure 7. 

1 . Magnetic Field Data 

Figure 8 shows the SCATHA magnetic field data 
recorded from 0300 UT to 0800 UT. Included are the total 
magnetic field, Bq, and the three topographic components of 
the field (Bjj,By,B 2 ). The significant rise in the magnitude 
of Bq, the total field, and B^, the northerly component, is 
due to the continuous descent of the spacecraft through the 
earth's magnetic field for that period of it's orbit. 
During that time period, SCATHA' s radial distance diminishes 
from 7.3 Re to 5.3 Re, and L decreases from 8.28 to 6.16. 
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NOON 



SCATHA DAY 210 ORBIT 




RE 



H FIRST OBSERVATION OF THE PCS 
e+ = MAXIMUM INCLINATION OF THE ORBIT 
r? H LAST OBSERVATION OF THE PCS 
€ s ORBIT CROSSES EQUATOR 
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e- H MINIMUM INCLINATION OF THE ORBIT 
s APOGEE 



Figure 7. Equatorial Plane Plot of SCATHA Orbit 
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DAY 210 MAGNETIC FIELDS 




Figure 8. Magnetic Field Data from 0300 to 0800 UT 
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GAMMA 



EAST-WEST MAGNETIC FIELD ( By ) 




Figure 9 . The Pc5 in the East-West Component 

of the Magnetic Field 
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The amplitude of the radial component of the 
magnetic field, varies from 1 to 4 gamma during this 
period. The magnitude of this component is largely due to 
the tilt of the earth's main magnetic field vector in the Z 
direction. The perturbations in are a result of a small 
component of the Pc5 in that direction. The Pc5 is observed 
in the plot of the east-west component of the field. By, 
beginning at 4:16 UT and continuing until just before dawn, 
06:56 UT. Since the perturbations in Bj^ and are small, 
the Pc5 event is largely contained in By. 

The By data are expanded in Figure 9 for the time 
period 0400 to 0700 UT (2:14 to 5:21 local time). A charac- 
teristic of this Pc5, which makes it fairly unique, is the 
packetized nature of the pulsations. These packets imply 
different regions of hydromagnetic resonance in the magneto- 
sphere. The packets are labelled A through E for con- 
venience in comparison to the electric field data which 
follows. The period of the pulsations decreases with time. 
It will be shown that the period is a function of L and of 
the particle density. 

2 . Electric Field Data 

Figure 10 is a plot of the electric field data 
measured by SCATHA from 0300 to 0800 UT. As mentioned 
earlier, the electric field data is despun into two com- 
ponents of the solar ecliptic (SE) coordinate system. The 
Ejj component points in the sunward direction and remains 
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DAY 210 ELECTRIC FIELDS 




Figure 10. Electric Field Data from 0300 to 0800 UT 
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almost constant at 1.5 mV/m throughout the event (2 mV/m has 
been added to the value of Ejj/ in Figure 10, in order to 
separate Ejj from Eg). This component is a result of the- 
SCATHA's relative motion through the plasma, the static 
convective electric field, the sunward electric field caused 
by photoemission [Ref. 17] and a partial projection of the 
Pc5 electric field perturbation vector in the X direction. 

The Eg component is perpendicular to the solar 
ecliptic plane. This component is nearly equal to the total 
electric field component of the Pc5. A -0.43 mV/m offset 
seen in the data is attributed to error fields and a steady 
state convective field. 

Figure 11a is an expanded plot of Eg from 0400 UT to 
0700 UT. The packetized nature seen in the magnetic field 
pulsation is reflected in the electric field. For con- 
venience, the electric field data is separated into named 
resonance packets as listed in Table 2. 

TABLE 2 

HYDROMAGNETIC RESONANCE PACKETS 

Packet t begin t end 

(hours UT) 



A 


4.41 


4.90 


B 


5.03 


5.46 


C 


5.65 


6.02 


D 


6.10 


6.40 


E 


6.63 


6.95 



The component of the electric field, Ep, which 
results from a perturbation in By (topographic coordinates). 
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Figure 11. Measured Electric Field and Low Energy Ion 
Distribution Function Showing Modulation 
of Ion Flux by the Field 



32 



should be perpendicular to both By and Bq. (SE 

coordinates) is believed to be a projection of Ep onto the 
electric field antenna. This is based on the assumption 
that there can be no component of the electric field 

sustained in the direction of Bq. 

6E • Bq = 0 (2.1) 

The electrons in a plasma are free to move along the 

magnetic field line and quickly reduce any potential 

difference which might be impressed along that direction. 

Figure 12 shows the geometry for extracting Ep from 
E 2 . Here a is the angle between Bj^ and Bq and also the 
angle between the radial direction and Ep. p is the sum of 
the satellite's geographic latitude and the satellite's 
inclination from the solar ecliptic plane. The angle $ 
between E 2 and Ep is given by 

$ = 90 - (a + p) (2.2) 

Ep - E 2 / cos $ (2.3) 

Table 3 lists these angles and the correction 
factors, e, for calculating Ep for the hours of 0400 to 0700 
UT. The worst case value is only 60 percent. The 
correction factor in Table 3, e, is used as follows: 



E = 1 / cos $ 
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a H ANGLE BETWEEN MAGNETIC FIELD LINE AND Bx 
<y = INCLINATION OF EARTH AXIS ON DAY BIO ( 1 B . 4 « ) 
0 = GtOGRAPHIC LATITUDE OF SCATHA SPACECRAFT 
^ " PLANE^^^^°'^ SATELLITE FROM SOLAR ECLIPTIC 

<3 = (e + tf”) at local MIDNIGHT 
© AT LOCAL DAWN 

$ = 90 - ( cx -V (3 ) 

Sp = Ez/cos $ 



Figure 12. Calculation of Ep 
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TABLE 3 

ELECTRIC FIELD DATA CORRECTION FACTORS 



(UT) 


a 


3 




e 






(degrees ) 






0400 


54.4 


19.2 


16.4 


1.04 


0430 


50.1 


18.0 


21.9 


1.08 


0500 


47.8 


16.7 


25.5 


1.11 


0530 


41.3 


15.2 


33.5 


1.20 


0600 


37.6 


13.4 


39.0 


1.29 


0630 


32.9 


11.5 


45.6 


1.43 


0700 


28.7 


9.2 


52.1 


1.63 



The factor which most affects the magnitude of this 
correction is the shape of the earth's magnetic field at 
different locations in the orbit. It can be seen in Figure 
12 that, for small a, less of the perpendicular field is 
projected onto the E 2 antenna. Near midnight the field is 
extremely elongated giving a large angle between Bq and B^, 
and thus a small angle between E and the Z antenna. At dawn 
the magnetic field most closely approaches a dipole and so 
this angle is small. During the analysis of the electric 
field, with this small correction in mind, the actual 
electric field data will be assumed to represent the perpen- 
dicular electric field (E 2 ~ Ep) • 

3 . Low Energy Ions 

Figure 11b shows a plot of the log (base 10) of the 
distribution function of the 39.7 ev ions measured from 4 to 
7 hours UT. The Pc5 can be seen in the ion flux. A 
comparison of this plot with the electric field shows that 
the 39.7 ev ion flux is modulated in phase with the 
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perturbation of the electric field, 6 E 2 . This reflects the 
convective drift induced in the ions by the electric field. 
Higher energy ions and electrons are also modulated, as 
shown in the expanded presentation of the D packet which 
follows . 

B. A CLOSER LOOK, THE D PACKET 

Since the D resonance packet exhibits the largest 
amplitude fields, it was chosen for closer analysis. This 
packet begins at 6.1 hours UT and ends at 6.4 hours UT. At 
this time SCATHA is at 7.3“ North, having recently passed 
the maximum inclination of the orbit. The D packet is 
observed for a period of 17.23 minutes. The center of this 
region is at L =6.97 and at 4:21 local time. During this 
time SCATHA 's average velocity was 3.25 km/sec. After 
taking into account the decrease in SCATHA 's radial position 
through this period, this resonance packet appears to span 
an arc of 3300 kilometers. 

1 . Field Relationships 

Figure 13 shows plots of spline fitted E 2 and By 
versus time for the D packet. Within the resonance region, 
E leads B by an average of 130® ±10®, (~ 3n/4), and the 
average period of the pulsations is 224 seconds. Averages 
are used since the periods and phases are continuously 
changing, as if the fields are driven by a multi- frequency 
or non-harmonic force. The electric field amplitude varies 
from 0.62 to 1.8 mV/m. The magnetic field amplitude varies 
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Figure 13. Pc5 Electric and Magnetic Fields 

(D Packet) 
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from 6 gamms to 12 gamma. Th© large amplitudes at the 
center of the packet distinguish the center of resonance in 
this region. 

The electric field vector and the magnetic field 
vector of the Pc5 lie in a plane perpendicular to the main 
magnetic field line. The electric field vector pulsates 
toward and away from the earth in a direction between the 
radial direction and the (SE) +Z direction. The magnetic 
field vector pulsates in a nearly east-west direction. If 
the fields were in phase the wave would be propagating in 
the -X direction (south along the magnetic field line). A 
standing wave will exhibit tt/ 2 phase difference. A phase 
difference of -tt would indicate the wave is propagating 
northward . 

The average phase of fields in Figure 13 is 3 tt/ 4. 
This implies that the wave is not a pure propagating wave. 
The continuously changing phase of the fields and the motion 
of th^ spacecraft through different regions of resonance 
make the phase interpretation somewhat ambiguous . 

The complexity of the phase relationship is born out 
by a phase hodogram plot of 6E2 versus 6By in the D packet. 
Figure 14. In a hodogram a propagating wave would give a 
diagonal line, and a standing wave a circle (if E and B are 
plotted in Volts/meter and Tesla). This hodogram implies 
that the wave is propagating at times, where the plot is 
nearly linear. At other times it shows the phase of a 
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Figure 14. Phase Plot of the Pc5 Electromagnetic Wave 

(D Packet) 
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standing wave, where the plot becomes more circular. The 
average phase of 3tt/4 is reflected in the overall elliptical 
nature of the plot. 

E and B are related by E = v^j^B. Hence, the slope of 
a line, or the inclination of an ellipse, in Figure 14 gives 
v<j). The slope of the major ellipse is 0.153 (mV/m) /( gamma) 
which gives a phase velocity of 153 km/sec. 

2 . Particle and Field Relationships 

Particle count rates were measured by the UCSD 
electrostatic analyzers and the energy distribution 
functions were calculated. Ion and electron fluxes are 
obtained for 63 energy levels every 16 seconds. 

Figure 15 contains plots of the log^^Q F versus time 
for low energy ions, high energy electrons, and two sets of 
high energy ions measured between 6.1 UT and 6.4 UT . The 
data sets chosen were from energy levels which exhibited the 
highest particle count rates. This was done in order to 
minimize the statistical error. 

A plot of log F versus time shows the fluctuations 
in the energy of the ions with time. The main cause of the 
change in the kinetic energy of the low energy ions is the 
EX B force, which is generated by the Pc5 wave. In the D 
packet region the magnetic field strength, Bq, is approxi- 
mately 140 gamma. Hence, a 2.6 mV/m electric field should 
give the plasma an E X B drift velocity equal to 

Vj 3 = 6 E 2 X Bq / |Bq 2 I ~ 18.6 km/sec (2.5) 
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Figure 15. Ion and Electron Distribution Functions 
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Here the upper limit of the electric field corrected with 
Table 3 values is used (6E • e). 

If the plasma flow is perturbed by a change in the 
fields, then the kinetic energy is given by 

U = 1/2 m (V - Vd)2 (2.6) 

= 1/2 m (V2 - 2 V Vj 3 cos 0 + Wjp-) 

where V is the particle and Vj^ the E X B flow. If V >> Vjj 
(which is true for 40 ev ions) and, if 0 = 0, then Equation 
(2.7) can be expanded. 

U = 1/2 m (V - Vd)2 (2.7) 

• 1/2 mV^ - (mV vjj) 

If the particle energy distributions were 
Maxwellian, they could be represented by 

F = n (m/27TkT)^/2 exp -(U/kT) (2.8) 

and by Equation (2.7) 

log F = constant + (log e) (mV/kT) • Vq 

Hence, the plot of log F for 40 ev ions is propor- 
tional to Vj). Comparision of 39.7 ev ions with the electric 
field in Figure 15a shows this to be the case. The ion flux 
is out of phase with the electric field by tt (-E 2 is 
plotted). The phase difference with the electric field is 
due to the unidirectional counting of the ESA's. Since they 
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are looking in the direction of the positive EX B force, 
the count rate is maximum when 6 E 2 is negative. Thus, the 
low energy ion flux is modulated by the change in the flow 
energy induced by the E X B force. The higher energy 
particles respond primarily to B [Ref. 22]. 

Figure 15b shows that 9.3 kev ions are effectively 
modulated in phase with the magnetic field data. These 
higher energy particles are, of course, unaffected by the 
small' change in the velocity caused by the E X B force. 
Hence, the changes in 9.3 kev ion flux reflect the effective 
motion of the different magnetic field lines with respect to 
the satellite. Figure 15c shows 4.7 kev electrons which, in 
like manner, show flux modulations in phase with the 
magnetic field. 

At higher energies the modulation of the ion flux 
begins to lag the magnetic field . This is shown in Figure 
15d where the flux of the 47.4 kev ions lags the magnetic 
field by an average of tt/ 2. This difference reflects the 
large gyroradius (225 km) of these high energy particles. 

These comparisons illustrate the point that the 
plasma is moving with the magnetic field lines of force. 
Hence, the plasma and the field lines oscillate together as 
if the particles are frozen to the field lines [Ref. 6 ]. 

3 . Plasma Parameters 

The plasma density, temperature and flow velocity 
can be determined by assuming that ion distribution 
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functions are nearly Maxwellian and fitting Maxwellian 
distributions to the data. 

Analysis of the ion data requires a multi-species 
population due to the various masses. The distribution 
function becomes 

F = S Fi = S n^ (mi/ 2 TTkT) 3/2 exp -(Ui/kTi) (2.9) 

and log F for each population is 

log F = log ( n (m/2TrkT)3/2 j _ (U/kT) log e (2.10) 
= f(p,T) - (log e / kT) U 
= f(p,T) + f(T) U 

To begin the fitting process the distribution 
function of the ion data is calculated and a plot of log F 
versus energy is made. As an example. Figure 16 shows a 
plot of the ion distribution function calculated from ion 
flux recorded by the N/S detector at 6:12 UT. This time was 
selected because it is near a large amplitude in the 
electric field. It is here that the highest drift 
velocities can be expected and the effect of the positive 
satellite potential shadowing the lowest energy ions (H"'’) 
the smallest. 

In this fitting process the plasma was assumed to 
contain H"*", He'*' and O"'' ions only. Once the ion species are 
assumed, there are three unknowns for each species in the 
distribution function. The temperature, kT, the density, n. 
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Figure 16. UCSD Ion Distribution Function 
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and the flow velocity, v. A self consistent set of values 
is required. In particular, consistent flow velocities are 
required since all ions move together. One additional free 
parameter is the satellite potential. 

The temperature, kT, is approximated to be 10 ev. 
In a more distinctly separated energy distribution, each ion 
species would show a well defined peak and the width of the 
peak would define the temperature. In the present case, the 
curvature of the distribution function is a result of both 
temperature and velocity. 

Since the detector cannot be assumed to be looking 
in a direction exactly parallel to the flow, a look angle, 
4), must be included. The choice of this angle affects the 
calculated density and flow velocity. In this case 10 
degrees was chosen. 

A peak in the distribution function determines a 
velocity, vp, given by 

Vp = (2U/m)^'^^ 

= vq cos <t> (2.11) 

where Vj^ is the plasma drift velocity. Hence, the error 
introduced in the calculated velocity by an erroneous 
selection of the view angle is proportional to cos <t>. An 
error of 15° in the approximation of <j), will cause an error 
in Vp of less than 10 percent. Olsen discusses this process 
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and shows that the error, introduced in the calculation 
of n is more severe [Ref. 23]. 

5 « exp cos^<t>) (2.12) 

In Figure 16, the peaks occur at approximately 8 ev, 
35 ev and 100 ev. The lower energy peak is assumed to be 
due to H"'', the middle energy peak. He'*', and the highest 
energy peak , O'*'. An approximation, Vp, for the drift 
velocity of each population can be made from the measured 
energy peak. 

Since V = (2U/m)l/2 (2.13) 

Vp 39 km/s 
Vp jjp • 41 km/ s 
Vp Q ~ 35 km/s 

In the example. Figure 16, the data are fitted with 
a Maxwellian distribution function for H'*’, for He'*' and for 
O’*". Consistency requires that the Maxwellian plots for all 
three species have the same drift velocity and a reasonable 
density ratio. The values used for this fit are 35 km/s 
velocity, a 10 ev temperature, a look angle equal to 10 
degrees, and ion number densities 1.0, 0.5, and 0.3 .(cm”^) 
for H'*', He'*' and O'*'. The results of this fitting process 
using UCSD data over the D packet are listed in Table 4. 
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TABLE 4 

PLASMA ION NUMBER DENSITY (1-500 ev ) 



Ion 


n 


V 




( cm" 3 ) 


(km/ sec ) 


H+ 


1,27 ± .67 


35 ± 5 


He+ 


0.38 ± .05 


35 ± 5 


0+ 


0.30 ± .10 


35 ± 5 


Total 


1.95 ± .73 


- 35 



The flow velocity of 35 km/sec contains both the flow 
induced by the PC5 and the steady background flow. It is 
estimated that the background convective flow is 1 - 10 
km/sec. The steady state electric field responsible for 
such a flow (0.2 - 2.0 mV/m) is poorly measured by the 
electric field instrument. The remaining, pulsation induced 
flow, is approximately 25 km/sec. 

From Equation (2.6), the 6E2 X Bq velocity was 
calculated to be 18.6 km/sec. This implies a steady state 
flow velocity near 16.4 km/sec. Since 6E = vjj • Bq and Bo = 
140 gamma, a convective electric field as large as 2.3 mV/m 
is implied from the electric field data. Analysis 
throughout the D packet is consistent with this interpre- 
tation. 

An alternative interpretation is that 5E is twice 
the value measured by the floating probe. Experience with 
the electric field data by the principal investigator, T. L. 
Aggson, has suggested problems with the amplitude of the 
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measured fields. An aggreement within a factor of two, at 
this point, is considered to be supportive of the direct 
measurements of E. 

The complete ion distribution is needed to determine 
the total ion density. The high energy UCSD data were 
combined with the high energy ion composition data from the 
SC8 High Energy Ion Composition Experiment [Ref. 21]. SC8 
determined average ion densities in the 0.1 to 32 kev range 

for the one hour period from 0600 to 0700 UT. The average 

density was (1.974 cm“^}, the He+ density was 
(0.089 cm“^) and the 0“^ density was (0.398 cm"^). The SC8 

data are shown in Figure 17. The line plot in Figure 17 is 
a plot of log F for the total population. The low energy 
(100 eV - 500 eV) O’*" ions are the same population inferred 
by the fitting process of the UCSD data in Figure 16. The 
UCSD (0 - 80 kev) ions are shown in Figure 18. The 10 keV 

ions are part of the bulk of the plasma sheet population, 

while the 47 keV ions are from the tail of the distribution 

function. The SC8 data are combined with the parameters 
obtained from the UCSD experiment and reported in Table 5. 

TABLE 5 

PLASMA ION DENSITIES (.01 - 32 KEV) 

Ion n 



( cm" ^ ) 



H+ 


3 . 24 


± 


.68 


He+ 


0.47 


± 


.15 


0+ 


0.40 


± 


.17 


TOTAL 


4 . 1 


± 


1. 
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Figure 18. UCSD Ion Distribution Function 
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4 . Electron Density 

The calculation of the ion density depends largely 
on the distribution function of the ion population being 
nearly Maxwellian. An error in the look angle can introduce 
a large error in n. Since the number density of the 
electrons should be equal to the number- density of the ions 
in a plasma, the value of ng is a good check of the ion 
density calculations. 

All the members of the electron population have the 
same mass, and that mass is known. Hence, one of the 
unknowns that hindered the ion calculation is removed . A 
useful characteristic of magnetospheric electrons is that 
sections of the distribution function plot have clearly 
defined slopes which identify different energy populations 
of electrons. 

In a plot of the electron distribution function 
versus energy 

log Fg = log {ng (mg/ 2 TrkTg) j _ igg g (i/kTg) Ug 

= c + b Ug (2.14) 

where b is the slope of each section of the distribution 
function and c is the y intercept. From Equation (2.14) 



kTg = -log e / h 


(2.15) 


ng = 10^ / (mg / 2 tt kTg)^/^ 
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Figure 19 is a plot of the electron distribution 
calculated from UCSD data taken at 6:12 UT. The units of 
Figure 17 are (sec^/km^) versus (keV) . Correcting for units 

(2.16) 

n© = (10^ X 10"24) / (9.05 x lO"^^ / (ev) )3/2 

In Figure 19 two energy sets of electrons are identified. 
By the slope intercept method, the 17.5 kev electrons are 
shown to have a number density of 0.007 cm“^. The number 
density of the 1.65 kev electrons is shown to be 1.77 cm“^. 
By expanding the lower energy section of this data two other 
electron groups were identified. The number densities which 
were calculated for each set of electrons is reported in 
Table 6 . 



TABLE 6 

ELECTRON NUMBER DENSITY 
Energy Range Temperature n 



(kev ) 


(kev ) 


( cm 2 


) 


15 - 80 


17.50 


0.007 


± 


.001 


1.0 - 15 


1.65 


1.77 


± 


.09 


0.2 - 1.0 


0.71 


1.46 


± 


.09 


0.0 - 0.2 


0.12 


0.67 


± 


.04 




TOTAL 


3.9 


± 


.2 



Hence, the calculated ion number density (4.1 cm“^) 
is consistent with the independently determined electron 
density. 
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Figure 19. UCSD Electron Distribution Function 
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5 . Self Consistency of the Data 

It is now important to determine whether these 
numbers represent a self-consistent relationship between n, 
B, E and V. Since the plasma acts as a parallel resistive 
circuit on the electric field antenna, a major concern is to 
verify the amplitude measurements of . 

Cummings used the following equation to deduce the 
density of the ambient plasma from field and particle data 
measured by the ATS-6 spacecraft in geosynchronous orbit 

n = 6 b 2 / <m> V£)2 (2.17) 

where <m> is the average particle mass [Ref. 7]. This 
relationship is used to investigate the consistency of the 
ion density calculations. 

Assuming that the plasma is made up of H"*", He"*" and 
O"*", and using the values from Table 6, the ion mass density 
is given by 

p = S nj^ mj^ (2.18) 

= (3.24(1) + .47(4) + .40(16)) mp 
p = (11.52 mp) kg/cm^ 

and <m> = p / n = 2.81 mp 

where mp = 1.672 x 10"^^ kg. 
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Equation (2.17) is then rewritten. 



= 5B^ / Pq <m> n 

V = 6B / (Po 

Since | v | = | 6 E / Bq | 

6E = 6B Bq / (Po p)l/2 (2.19) 

= 6.4 X 10"3 (5B Bq) mV/m 
where the units of 6B and Bq are gamma. 

In the D packet Bq = 140 gamma 

6 < 6By < 12 gamma (2.20) 

Therefore the magnitude of the electic field inferred by 
mass density is 

5.4 < 6E <10.8 mV/m (2.21) 



The measured values of 6E2 were 

0.62 ^ 6E2 ^ 1.8 mV/m (2.22) 

Taking into account the electric field correction factor of 
1.43, from Table 4, the magnitude of the perturbed 
electric field is calculated to be 

0.9 < 6Ep < 2.6 mV/m (2.23) 

These values are a factor of 3 lower than predicted by 
Equation (2.21). A factor of 2 could be accounted for if 
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the impedance of the plasma was comparable to the internal 
impedance of the electric field experiment (this is roughly 
the case). A small error, not larger than a factor of 2, in 
the electric field data is implied by the comparison. 

6 . The Alfven Velocity 

The same comparison can be made in terms of 
velocities. The Alfven velocity is 

= co/k = Bo / (Po (1-2) 

Curl E = -dB/dt 
ik Ep - -(-ico) By 

6Ep/By - w/k = (2.24) 

Using upper and lower bounds of the field amplitudes 
given by Equations (2.20) and (2.23), the Alfven velocity 
calculated from measured E is 

328 km/sec < < 378 km/sec (2.25) 

Using the ion density relationship 

Va = Bo / (Mo = 900 km/sec (2.26) 

The electric field data appears to be low by a factor of 3 
in this comparison. 

7 . The Poynting Vector 

The nature of the propagation of the PCS can be seen 
by examining the relationship of 6 E 2 and 6By. If the 



Since 

then 

and 
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micropulsation was a classical plane Alfven wave, then 

S = (1/y) E X B 

- (l/Po) 5E X 6B 

- 1 (l/Wo) (SEy 6 B 2 - 6E2 6 By) 

which assumes the other components of the Poynting vector 
are small. 6 Ey is not measured, but 6 B 2 is small, so 

S - 1 (l/Po) (-5Ez 6By) (2.27) 

In Figure 20 Px = ( 5 E 2 6By), so 

S = -1 (1/jiQ) (Px) X 10“^^ Watts/m^ (2.28) 

The negative values of Px in Figure 20 indicate that 
the wave is propagating in the +X direction, northward along 
the magnetic field line. The wave must have been generated 
south of 7° N. The fact that Px is, at times, zero or 
positive indicates that the micropulsation is not a simple 
traveling wave, but is superimposed on a smaller wave 
traveling in the southward direction. This is the mixed 
standing / propagating nature inferred previously in Figure 
14. 

And finally, the magnitude of the Poynting vector is 
calculated 

Px = (1.8 mV/m) (-12 gamma) 

|S| = (l/Po) |Px| X 10“12 (2.29) 

= 17.2 ywatts/m2 
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This presumes that the electric field measurement is 
correct. Given the possible error in 6E, the magnitude of 
the Poynting vector could be as high as 50 yWatts/m^. 

8 . Dependence of the Period on L 

The period of the Pc5 depends on the length of the 
field line and the Alfven velocity of the wave [Ref. 24]. 
The variation of micropulsation period with latitude can be 
expressed as a function of L and n, the plasma number 
density. [Ref. 25]. 

•c = (2tt Re/Be) * (Po n/Q)l/2 (2.30) 

Bg = 3.12 X 10“^ Tesla 

Q is proportional to an Alfven wave eigenvalue which is 
calculated for different values of L, density gradient, and 
wave type. In a simpler form this equation corresponds to 

T = 2tt (m/k)l/2 ^ 

The SCATHA data can be used to determine the dependence of 
the period of this Pc5 on L. If the period, x, is 
proportional to 

log X = log C + m log L (2.31) 

where C is a constant. 

Figure 21 shows the average micropulsation period in 
each resonance packet, plotted against L. The slope, 
determined by linear regression, is 3.96 ± .07. Hence, the 
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Figure 21. Pc5 Period Versus Mcllwain L 
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SCATHA data is in close agreement with the theoretical 
relationship 

X “ (2.32) 

Having verified the L dependence of the period. Equation 
(2.30) is rewritten and utilized to predict a mass density 
to compare with the previous measurements . 

p = n <m> = (X Bg / 2tt Re )^ X (Q / Pq) (2.33) 

For the values x = 224 sec (in the D packet) 

L = 7.0 

Be = 3.12 X 10“^ Tesla 

n <m> = 2.514 mp • 0 (cm”3) (2.34) 

Orr and Matthew calculate Q for a toroidal mode wave 
in plasma density gradients which obey power laws between 
and r 2 [Ref. 25]. At L = 7 these values are 

3.6 (for R'^) < Q < 4.2 (for R^) (2.35) 

The mass density inferred by Equation (2.33) is 

9,05 mp < n <m> < 10.56 mp (cm“^) (2.36) 

The measured mass density was 11.52 mp (cm“^). If the 
theoretical relationship is valid for this Pc5 (it has been 
shown that it is substantively valid by the agreement with 
the power law of the period), then the difference between 
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predicted particle density and measured density implies that 
the mass density has been slightly overestimated. However, 
a deviation of the wave period of only 10 seconds from the 
average period, used in Equation (2.34), would bring the 
predicted mass density into exact agreement with the 
measured mass density. 

Since the values for O were derived from a 
theoretical eigenvalue equation for the toroidal wave, 
complete aggreement with the mass density results of this 
study is not demanded. However, this comparison shows 
strong agreement with both variables, L and p. The plasma 
density gradient is the unknown variable. 

C. SOLAR AND MAGNETIC FIELD CONDITIONS 

Geomagnetic and Solar Data, published by the Journal of 
Geophysical Research, indicates that Day 210 was a 
moderately active day, Kp = 3 to 5, with no magnetic 
substorms recorded. Forty eight hours prior to the Pc5 
event there was a period of very high magnetic activity in 
which a 12 hour magnetic storm occurred [Ref. 26]. These 
conditions support the statement by Jacobs [Ref. 1] that 
Pc5 ' s generally occur during recovery periods following 
magnetic disturbances. 

D. OTHER OBSERVERS 

SCATHA observed the Pc5 at 331 degrees east longitude on 
L = 7 . Magnetograms from the National Institute of Polar 
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Research from the conjugate ground stations at Leirvogur, 
Iceland (68.2* N 338.3* E) and at Syowa Base, Antarctica 
(69.0* S 39.6* E) indicate that this event was observed at 
these stations simultaneously from 0600 to 0700 UT. The 
time resolution of the ground station data was not clear 
enough to say more. These stations are located at the 
invariant latitude for L = 7 . The simultaneous observation 
of this event would imply that the generation was occuring 
near the equator. 

There is no ' indication of this observation at Narssarss- 
uaq Station, Greenland (61.2* N 314.2* E) or at St. Johns 
Station (47.6* N 307.3* E) . Magnetograms from the NOAA 
GOES 2 and GOES 3 satellites were also scanned. They showed 



frequent magnetic 


activity , 


but 


none 


which could 


be 


correlated with 


the 


Pc5 . 


These satellites were 


in 


geosynchronous orbit 


at 


252* 


East 


and 225* 


East 



respectively . 

The lack of Pc5 activity at locations west of the SCATHA 
sector indicates a limit on the temporal and spatial extent 
of the event. The Pc5 appears to have been confined to 45 
degrees of longitude in the morning sector observed by 
SCATHA. 
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III. DISCUSSION 



A. THE CHARACTER OF THE PCS 

In Chapter II, the observable quantities characterizing 
this Pc5 event were investigated by focusing on data taken 
at the temporal and spatial center of the event. The wave 
was shown to be a traveling Alfven wave guided by the 
magnetic field line. The Alfven velocity was nearly 900 
kilometers per second. A calculation of the Poynting vector 
showed the wave was traveling northward. 

The electromagnetic field perturbation was shown to be 
occuring in a plane perpendicular the the magnetic field 
line. The magnetic field lagged the electric field by 3 tt/ 4. 
The amplitude of the magnetic field perturbations was near 
12 gamma. The measured electric field reached 2.57 
millivolts per meter. By requiring consistency of the data, 
the measured electric field was shown to be low by a factor 
of 2 to 3 . The period of the Pc5, at L=7, was 224 seconds 
and the period of the wave was shown to be proportional to 
l4. 

The plasma density was estimated to be 1.9 x 10“^^ 
(gm/cm^). The ion number density was determined to be 4.1 
(cm“2). This is consistent with previous magnetospheric 
studies. The calculation of the Alfven velocity ~ 900 
km/sec is also consistent with previous studies. Most 
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recently, empirical models of average magnetospheric plasma 
and field data were developed to study Alfven and Fast Mode 
wave speeds [Ref. 27]. These models predict plasma 
densities of 3-5 (cm~^), plasma temperatures of 30-100 ev 
and Alfven wave speeds of 500 - 1000 km/sec at L = 7.0 in 
the morning sector of the magnetosphere. 

B. SPATIAL EXTENT OF THE PC5 

The most interesting factor in this Pc5 event is the 
packetized nature of the micropulsation resonances through 
this sector of the magnetosphere. The breadth and location 
in these regions is presented in Table 7 . 



TABLE 7 

RESONANCE PACKET LOCATION AND SIZE 



t 

(UT) 


ID 


Width 

(km) 


Latitude 
(“ N) 


L 


Local 

Time 


r 


4.66 


A 


5300 


7.38 


7.8 


2:45 


0.45 


5.25 


B 


5300 


7.50 


7.5 


3:23 


0.50 


5.84 


C 


4200 


7.45 


7.2 


3 :54 


0 . 32 


6.26 


D 


3300 


7.29 


6.9 


4:21 


1.00 


6.80 


E 


3600 


6.96 


6.7 


4:55 


0.72 


7.46 


F 


5400 


6.11 


6.4 


5:42 


0.49 



The r factor, in Table 7, is a normalized resonance 
strength for each region, which is the ratio of the electric 
field amplitude in a region to the electric field in the D 
packet. These resonance regions are seperated by sectors of 
low resonance which vary in distance from 1200 kilometers to 
3600 kilometers. 
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C. UNRESOLVED ISSUES 



A complete understanding of this event would require a 
vector description of the data. There are significant 
aspects of the Pc5 which were neglected by this study. One 
is the effect that the nutation of the spacecraft spin axis 
has on the raw data. The nutation effect can be seen as a 
small pulsation in each component of the magnetic field data 
in Figure 8. Another aspect is the significance of the 
smaller component of the Pc5 which is in the radial 
direction, The component has an average amplitude of 
4 gamma. This component has a small effect on the amplitude 
of the Pc5 magnetic field, but has a significant effect on 
the rotational motion of the field . A study that included 
this component would provide a more realistic description of 
the Pc5 event. 

A most difficult aspect, which is left for the student 
who follows this study, is the separation of the spatial 
character of the event from the temporal description which 
this study provides. Many of the effects seen in the data, 
such as the sinusoidal variation of the resonance in each 
packet, have been described but not explained. 

D . RECOMMENDAT I ONS 

1 . The spatial and temporal extent of this Pc5 could 
not be determined by SCATHA alone. Specifically, it is not 
known if local dawn was a boundary of the Pc5 or that the 
Pc5 ended at all locations at 6:56 UT. The radial boundaries 
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of the event are also left unresolved. The observation of 
the Pc5 ended when SCATHA passed L = 6.6. This may have 
been a coincidence, but it could be a major characteristic 
of this phenomenon. A correlation with a satellite that was 
East of SCATHA could provide an answer (such data may be 
available from the European GEOS-2 satellite). Better time 
resolution of the conjugate ground station data would 
determine the location of the generating mechanism. 

2. Approximate densities were calculated in the region 
of highest resonance, the D packet. The wealth of 
information made available by SCATHA provides an area for 
further study into the magnetospheric plasma profile between 
5 . 5 and 7.5 Re . 

3. The theoretical models presently available describe 
standing and compressional Alfven waves. Since this event 
provides extensive data from a traveling wave event, further 
investigation into this area of micropulsation theory is 
warranted . 

4 . Theorists have recently shown great interest in 
understanding the resonance conditions for Alfven waves. As 
emphasized in this paper, the packetized nature of the 
resonance conditions in the morning sector of the magneto- 
sphere is displayed vividly by the SCATHA data. These 
observations should be used to improve the understanding of 
Alfven wave resonance. 
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IV. 



CONCLUSIONS 



Observations of a Pc5 event occuring between L = 6.6 and 
L = 8.0 in the pre-dawn sector of July 29, 1979, at 7® N, 
were reported. Hydromagnetic resonance regions spanning 
3000 to 5000 kilometers with centers separated by 2400 
kilometers were identified. 

The plasma density in the highest resonance region was 
shown to be 1.92 x lO”^^ gm/cm^ . ion flux at 40 ev was 
shown to be in phase with the electric field . Ion fluxes 
at 9 kev, and 5 kev electrons, were shown to be in phase 
with the magnetic field. The 47 kev ions lagged the 
magnetic field by tt/2 radians. 

Particle data from the UCSD electrostatic analyzers was 
shown to be an effective means of cross checking electric 
field data. The electric field data was assumed to be a 
projection of the perturbation electric field onto SCATHA's 
antenna. By comparing this data with ion density 
predictions, the amplitude of the measured electric field 
data was shown to be a factor of 2-3 too low to be 
consistent with the complete set of data (this is considered 
to be a small error) . 

The periods of the micropulsations, between L = 6.3 and 
7.7, were shown to be proportional to the fourth power of L, 
which supports present models of resonant Alfven wave 
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activity in the magnetosphere. Alfven velocities were shown 
to be nearly 900 kilometers per second within high resonance 
regions of the event. 

The electromagnetic pulsations were in a plane 
perpendicular to the magnetic field line. In the Pc5, the 
magnetic field lagged the electric field by 3 it/ 4. The 
Poynting vector was shown to be in the direction of the main 
magnetic field vector, confirming that the wave was a 
traveling Alfven wave. Simultaneous observations by 
conjugate ground stations indicated that the generation 
could be occurring at the equator. 

The unique packetized resonance regions exhibited by 
this PCS as well as the non-classical wave character of this 
event indicates that further theoretical study utilizing 
data from SCATHA is warranted. 
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